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Abstract. The present understanding of nuclear electromagnetic properties including 
electromagnetic moments, form factors and transitions in nuclei with A < 10 is reviewed. 
Emphasis is on calculations based on nuclear Elamiltonians that include two- and three- 
nucleon realistic potentials, along with one- and two-body electromagnetic currents de¬ 
rived from a chiral effective field theory with pions and nucleons. 


1 Introduction 

A major goal in nuclear physics is to understand nuclear structure and dynamics in terms of under¬ 
lying interactions occurring between individual nucleons. Studies grounded on this basic picture of 
the nucleus are referred to as ab initio. An exceptionally powerful tool to asses the validity of our 
theoretical models is to investigate nuclear electromagnetic (e.m.) observables, such as ground state 
properties, e.g., e.m. moments and form factors, as well as e.m. reactions, e.g., photo- and electro- 
induced reactions. In these processes, external e.m. probes interact with the nuclear charge and 
current distributions with a strength characterized by the fine-structure constant a ~ 1/137. The small 
value of the fine-structure constant allows for a perturbative treatment of the e.m. interaction, while 
non-perturbative physics pertain only to the nuclear target. For light nuclei, terms that go beyond the 
leading order contribution in the Zcr-expansion (where Z is the number of protons) can be safely dis¬ 
regarded, leaving us with relatively simple reaction mechanisms and manageable formal expressions. 
For example, at leading order, the cross section associated with inclusive electro-nucleus scattering 
processes is factorized into the leptonic tensor, which is completely specified by the measured elec¬ 
tron’s kinematic variables, and the hadronic one associated with the nuclear target, and proportional to 
matrix elements squared of the nuclear e.m. charge and current operators. A clear connection between 
measured quantities, i.e., cross sections, and calculated matrix elements is then realized. Experimen¬ 
tal data of e.m. observables are, in most cases, known with great accuracy providing us with viable 
and strong constraints on our models. Likewise, for light nuclei, theoretical calculations are affected 
by relatively small statistical errors because for these systems the many-body problem can be solved 
exactly or within controlled approximations. This allows for solid comparisons between experimental 
data and theoretical predictions. 

In Fig. [7] a cartoon picture of the double differential cross section for electron scattering off nuclei 
is represented. Different values of energy u> transferred to the system, correspond to different excita¬ 
tion energies of the nucleus. By varying co, we can access the ground state (elastic peak), low-lying 
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Figure 1 . Schematic representation of the double differential cross section for electron scattering off nuclei at 
fixed value of momentum transfer. 

(discrete) nuclear excited states, giant resonance modes, and the quasi-elastic energy region up to the 
pion-production threshold. For each value of excitation energy u>, one can study the matrix elements’ 
behavior as a function of the momentum |q| transferred to the nucleus. In particular, by varying |q| 
one can explore the e.m. charge and current distributions with a spatial resolution oc 1 /|q|. In this talk, 

1 will focus on ab initio calculations of ground state nuclear e.m. properties, that is e.m. moments 
and elastic form factors, as well as widths of e.m. transitions occurring between low-lying nuclear 
states. These studies have been recently reported in a topical review on e.m. reactions on light nu¬ 
clei 01, where more details and references to original articles can be found. Recent developments 
on theoretical ab initio investigations on other very interesting e.m. processes in light nuclei, such 
as photo-absorption and radiative capture reactions, Compton scattering, sum rules ..., are well rep¬ 
resented in this conference, see, e.g., contributions by X. Zhang, J. Dohet-Eraly, H. Griesshammer, 
M. Miorelli, S. Bacca, N. Barnea, D. Rozpedzik, and A. Lovato in these proceedings. 

A theoretical understanding and control of nuclear e.m. structure and dynamics is a necessary 
prerequisite for studies on weak induced reactions, such as neutrino-nucleus interactions. The ex¬ 
perimental data acquisition for this kind of processes is comparatively more involved owing to the 
tinier cross sections and to the fact that neutrinos are chargeless particles and, thus, they are hard to 
collimate and detect. An important advance in this direction has recently been carried out by Lovato 
and collaborators 0], and for a status report on ab initio calculations of weak response functions in 
4 He and 12 C I refer to the plenary talk of A. Lovato (the associated contribution can be found in these 
proceedings). Moreover, a theoretical understanding of the structure and dynamics of light nuclei 
is a necessary prerequisite for research projects aimed at studying larger nuclear systems. For these 
reasons, it is imperative to first validate our theoretical understanding of e.m. reactions on light nuclei. 

2 Nuclear Hamiltonians and electromagnetic currents 

In the ab initio framework, the nucleus is described as a system made of A non-relativistic point-like 
nucleons interacting among each other via many-body forces and its energy is approximated by the 
following Hamiltonian: 



( 1 ) 


where K) is the non-relativistic single-nucleon kinetic energy, while v,j and F, ; / f are two-nucleon (NN) 
and three-nucleon (3N) potentials, respectively. Implicit in the equation above is the assumption the 
four-nucleon forces and higher order terms in the many-body expansion are suppressed. The NN 
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Figure 2. Diagrams illustrating one- and two-body EM currents entering at LO {e QT 2 ), NLO (e Q -1 ), N2LO 
(e Q°), and N3LO (e Q *). Nucleons, pions, and photons are denoted by solid, dashed, and wavy lines, respec¬ 
tively. 


and 3N potentials are phenomenological in nature in that they involve a number of parameters— 
subsuming underlying Quantum Chromodynamics (QCD) effects—that are fixed by fitting experi¬ 
mental data. For example, NN potentials are constrained to reproduce a large number of NN scattering 
data, along with the deuteron binding energy. Nuclear forces belonging to this class of highly accu¬ 
rate nuclear potentials are referred to as ‘realistic’. Most realistic potentials describe the long range 
(oc l/m„ where m„ is the pion mass) part of the nuclear interaction in terms of one-pion-exchange 
interaction mechanisms. Different dynamical schemes are implemented to account for intermediate 
and short range effects, among which multiple-pion-exchange, contact interactions, heavy-meson- 
exchange, or excitations of nucleons into virtual A-isobars. Here, the realistic potentials utilized to 
solve the Schrodinger equation //|'F) = E |'F) (where V F) is the nuclear wave function) are the Ar- 
gonne «18 0 (AV18) NN potential in combination with either the Urbana IX 0] or Illinois-7 0 3N 
potentials, as well as combinations of NN and 3N potentials derived from chiral effective field theory 
CfEFT)®. 

Nuclear charge ( p ) and current (j) operators describe the interactions of nuclei with external e.m. 
probes. They are also expanded in a series of many-body operators as 

P = 2jPi( q ) + ^Pi/q) + --- . (2) 

i i<j 

j = 2jji( q ) + ^j»v( q ) + --- . 

i i<j 

where q is the momentum transferred to the nucleus. In Impulse Approximation (IA), that is retaining 
only one-body operators in the equations above, nuclear e.m. charge and current distributions are 
simply the sums of those associated with individual protons and neutrons. The non-relativistic charge 
operator for point-like nucleons is simply the proton charge, while the nucleon current consists of a 
convection term associated with the current generated by moving protons and a spin-magnetization 
term associated with the spins of both protons and neutrons. The IA picture of the nucleus is, however, 
incomplete as it fails to explain, e.g., the measured magnetic moments of light nuclei. Corrections 
that account for processes in which external e.m. probes couple to pairs of interacting nucleons, de¬ 
scribed by two-body current operators, need to be incorporated in the theoretical ab initio description. 
Meson-exchange currents (MEC)—postulated in the ’40s by Villars id and Miyazawa d — follow 
naturally once meson-exchange mechanism are invoked to describe interactions between individual 
nucleons. They account for processes in which the external e.m. probe couples with mesons being 
exchanged between nucleons. The first evidence of meson-exchange effects in light nuclei can be 
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traced back to the 1972 work by Riska and Brown [12], in which MEC were found to provide the 
missing 10% correction to the IA value necessary to reach agreement between the calculated and the 
measured cross sections for the radiative capture of proton on neutron at thermal neutron energies. 
Since then, MEC have evolved into highly sophisticated and accurate currents. In their most recent 
formulation min, in order to assure consistency between nuclear forces and e.m. currents, MEC 
are constructed from realistic NN and 3N potentials so as to satisfy the continuity equation. Addi¬ 
tion of these MEC corrections to the IA picture successfully explains a wide number of e.m. nuclear 
observables in light nuclei j ]_5. 16]. 

Recent years have witnessed the tremendous development and success of ^ EFT ddl that rein¬ 
forces and grounds the achievements of conventional theoretical approaches. The relevant degrees of 
freedom of nuclear physics are bounds states of QCD, i.e ., pions, nucleons, and A’s,.... On this basis, 
their dynamics is completely determined by that associated with the underlying degrees of freedom 
of quarks and gluons, that is QCD. However, at low energies, QCD does not have a simple solution 
because the strong coupling constant becomes too large and perturbative techniques cannot be applied 
to solve it. ^EFT is a low-energy approximation of QCD valid in the energy regime where the typical 
momenta involved, generically indicated by Q, are such that Q <sc A^ ~ 1 GeV, where A x is the 
chiral-symmetry breaking scale. ^EFT provides us with effective Lagrangians describing the interac¬ 
tions between pions, nucleons, and A’s that preserve all the symmetries, in particular chiral symmetry, 
exhibited by the underlying theory of QCD at low-energy. These effective interactions, and the tran¬ 
sition amplitudes derived from them, can be expanded in powers of the small expansion parameter 
Ql A x , restoring, in practice, the possibility of applying perturbative techniques also in the low-energy 
regime. The unknown coefficients of this expansion in small momenta—referred to as low-energy 
constants (LECs)—while being tied to QCD effects and therefore attainable from QCD calculations, 
are, in practice, fixed by comparison with the experimental data. Due to the chiral expansion, it is 
then possible to evaluate nuclear observables to any degree v of desired accuracy, with an associated 
theoretical error roughly given by (Q/ A a ,) (v+1 \ This calculational scheme has been widely utilized 
to study both nuclear forces and nuclear electroweak currents. The many-body operators emerging 
from direct evaluations of the transitions amplitudes with interactions provided by ^fEFT Lagrangians 
involve multiple-pion exchange operators, as well as contact-like interaction terms. Nuclear two- and 
three-body interactions were first investigated in the late ‘90s by Ordonez, Ray, and van Kolck using 
a^EFT with pions and nucleons lE(il-22l. Currently, chiral NN (3N) potentials commonly used in 
ab initio calculations include up to next-to-next-to-next-to leading order or N3LO (next-to-next-to 
leading order or N2LO) corrections in the chiral expansion m m. 

Vector e.m. currents_ have been first derived from a ^EFT with pions and nucleons by Park, 
Min, and Rho in Ref. 1241. The resulting operators account for two-pion exchange terms entering 
at N3LO in the chiral expansion. These currents have been utilized in a number of so called ‘hybrid’ 
calculation^ of nuclear e.m. observables, including magnetic moments and Ml properties of A = 2- 
3 nuclei and radiative capture cross sections in A - 2-4 systems 25142711. More recently, ^EFT 
e.m. currents and charge have been derived up to one loop contributions included within two different 


imp lementations of time ordered perturbation theory: one is by the JLab-Pisa group (see Refs. [28 
31IB and the other one is by the Bochum-Bonn group (see Refs. 1321. 13311 ). In this talk, I focus on results 
obtained using chiral e.m. currents, and compare, where possible, different theoretical evaluations 
against the experimental data. For results based on conventional e.m. MEC currents I refer the reader 
to the review articles of Refs. EHH and references therein. 


1 In this kind of calculations, matrix elements of the chiral e.m. charge and current operators are evaluated in between wave 
functions obtained from conventional realistic potentials, as opposed to potentials derived consistently from ^EFT. 
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Figure 3. Left: Deuteron charge form factor from Ref. 13 ill (magenta hatched band) and Ref. |[36] (purple solid 
bands) compared with experimental data (black points). Middle: Same as left panel but for deuteron quadrupole 
form factor. Right: Deuteron magnetic form factor from Ref. frill (magenta hatched band) and Ref. |37] (purple 
solid band) compared with experimental data (open circles). 

Before I proceed presenting applications to e.m. observables, I will briefly describe the e.m. 
operators as they emerge from a^EFT with pion end nucleons. I will start off with the vector e.m. 
current which is diagrammatically represented in Fig. [2] The leading order (LO) contribution to the 
e.m. current illustrated in panel (a) is simply given by the non-relativistic one-body current used in 
IA calculations, while the N2LO one-body operator of panel (d) is a relativistic correction to the LO 
IA current. Currents of one- and two-pion range, describing long and intermediate range dynamics, 
enter at NLO and N3LO—panels (b), (c), (e)—(j). Short-range dynamics is encoded by the contact 
currents of panel (k). Unknown LECs enter the tree-level diagram of panel (e) and contact currents 
of panel (k). LECs entering the contact terms are of two kinds, namely minimal and non-minimal. 
The former enter also the NN chiral potential at NLO, and can then be constrained to NN scattering 
data; the latter need to be fixed from e.m. experimental data. A common procedure implemented to 
reduce the number of unknown non-minimal LECs (there are 5 of them) is to impose that the two 
LECs entering the isovector part of the tree-level current illustrated in panel (e) are in fact saturated 
by the A-couplings entering the A transition e.m. current |29l 311. The remaining three LECs are 
commonly fixed so as to reproduce the magnetic moments of the deuteron, triton, and 3 He 131 ]_. 

Early investigations on the e.m. charge operator in^EFT have been carried out in Refs. 134 43611 . 
and more recently loop corrections have been derived in Ref. [32], and subsequently Ref. [30]. In 
closing, we note that the structure of the charge operator is quite different from that of the vector e.m. 
current. Two-body corrections, in this case, are expected to be relatively small. In fact, leading two- 
body operators of one-pion range are suppressed as they enter at N3LO (as opposed to NLO as seen 
in the case of the vector currents), while there are no free LECs entering the charge operator [ 


3 Deuteron, 3 He and 3 H electromagnetic form factors 


For A = 2-3 nuclei, theoretical calculations performed by different groups are available, which makes 
it possible to compare them not only with the experimental data but also between themselves to test 
the solidity of the ab initio prescription. The left and middle panels of Fig. [3j show the deuteron 
charge and quadrupole form factors, respectively, calculated by Piarulli and collaborators 13 111 (ma¬ 
genta hatched bands) and Phillips II35L13611 (purple bands). Both calculations are based on chiral NN 


potentials. In particular, Piarulli et al. use wave functions from the chiral NN potential at N3LO [38], 
while Phillips those from the NN interaction at N2LO [|39!]. The thickness of the bands represents 
the sensitivity of the results to different cutoffs utilized to regularize the divergent behavior at high 
momenta of the chiral operators’ matrix elements 131.35. 3(Jl■ The two calculations very nicely agree 


with the experimental data for low-value of momentum transferred (q ^ 3 fm 1 ) and exhibit a similar 
(small) cutoff dependence. In the case of the charge form factor, as q increases, both theoretical cal- 
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Figure 4. The 3 He and 3 H magnetic form factors, and their isoscalar and isovector combinations from Ref. 13 in 


culations exhibit a more pronounced cutoff dependence, and differ between each other, an indication 
that this observable is sensitive to the nuclear wave functions utilized in the calculations. In the case 
of the quadrupole form factor, the agreement with the experimental data is seen up to q - 6 fm 1 , 
well beyond the expected regime of validity of the^EFT framework. In the right panel of Fig. [3] we 
compare the results for the deuteron magnetic form factor obtained by Piarulli et al. 1311 (hatched 
magenta band) based on chiral N3LO potential [38] and chiral e.m. currents at N3LO 11291] . with the 
fully consistent ^EFT calculations by Rolling et al. [37] (solid purple band) based on the chiral NN 
potential at N2LO [39] and chiral e.m. currents at N3LO 132. 331. The theoretical results are in very 
good agreement with each other and with the experimental data for values of momentum transferred 
q 3 fm 1 , and present a comparable cutoff dependence. 


Form factor calculations of A — 3 nuclei have been reported in Ref. 13111 . Here, we show results for 
the trinucleon magnetic form factors obtained utilizing the chiral e.m. currents at N3FO of Refs. 128 . 
[29fl and two sets of nuclear Hamiltonians, namely the AVI8 fl] NN plus UIX 0] 3N potentials, and 
the N3LO 0 NN and N2LO § 3N potentials. Calculations in IA are given in light blue (based on 
chiral interactions) and blue (based on conventional interactions), while full calculations that include 
the complete e.m. current up to N3LO are given in magenta (based on chiral interactions) and red 
(conventional interactions). In the figure, the top panels show the 3 He and 3 H magnetic form factors, 
while the bottom ones show their isoscalar (F£) and isovector ( F x '.) combinations B31I1 . As it is well 
know from studies based on the conventional approach (see Ref. il3l), two-body e.m. currents are 
crucial to improve the agreement between the observed positions of the zeros and the predicted ones 
at FO (or IA). Despite the excellent agreement between theory and experiment for q < 2 far 1 , the 
theory underpredicts the data at higher momentum transfers, while the zeros are found at lower values 
of q than observed. The theoretical description of the first diffraction region is still incomplete. 


4 Magnetic moments and electromagnetic transitions in A < 10 nuclei 


Moving on to larger nuclear systems, we find a number of Greens’s Function Monte Carlo calcula¬ 
tions 1 41 1 based on the AVJ_8 jj] plus IL7 0 nuclear Hamiltonian that use the chiral e.m. currents up 
to N3LO from Refs. 128. 29. 3J |. Magnetic moments of light nuclei [40;] are summarized in the left 
panel of Fig. 0 where IA results are given by blue dots, while calculations that include the full chiral 
e.m. current operator are indicated by red diamonds to be compared with the experimental data rep¬ 
resented by black stars. First, we note that corrections from two-body currents are found to be small 
where the IA picture is satisfactory (see, e.g., 6 Li, 9 Be, 10 B), and large where the IA picture is incom¬ 
plete (see, e.g., 7 Li, 7 Be, 9 C, 9 Li). Corrections from two-body components can be as large as 40%, 
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Figure 5. Left: Magnetic moments in nuclear magnetons for A < 10 nuclei from Ref. Qjp]. Black stars indicate 
the experimental values, while blue dots (red diamonds) represent GFMC calculations which include the LO or 
IA (N3LO) e.m. currents from^EFT. Right: Transition widths from Ref. El normalized to the experimental 
values for A = 6-9 nuclei. The notation is as in left panel. 


as seen in the case of 9 C, and are crucial to reach (or improve) the agreement with the experimen¬ 
tal data. It is also interesting to note that two-body effects, while being significant for the 9 C and 9 Li 
magnetic moments, are found to be negligible for those of 9 Be and 9 B. This behavior can be explained 
considering the dominant spatial symmetries of the nuclear wave functions for these A - 9 systems. 
For example, the dominant spatial symmetry of 9 Be ( 9 B) corresponds to an [a, a, n(p)] structure [42]. 
Therefore, the unpaired nucleon outside the a clusters does not interact with other nucleons. As a 
consequence, two-body currents, that describe the coupling of external e.m. probes with pair of inter¬ 
acting nucleons, produce a negligible correction. On the other hand, the dominant spatial symmetry 
of 9 C ( 9 Li) corresponds to an [a, 3 He ( 3 H), pp {mi)] structure, and two-body correlations contribute 
in both the trinucleon clusters and in between the trinucleon clusters and the valence pp {mi) pair, 
resulting in a large two-body current contribution. 


GFMC calculations of selected E2 and Ml transitions in low-lying nuclear states 0 are summa¬ 
rized in the right panel of Fig.0 Predictions in IA are represented by blue dots, while those obtained 
with the full chiral e.m. current operator are represented by red diamonds. Calculations for E2 tran¬ 
sitions implicitly include the effect of two-body currents via the Siegert theorem, where the charge 
density is used in IA. Also for these observables the effect of two-body e.m. currents can be large, and 
for cases in which the experimental errors are relatively small, e.g., 7 Li(l/2 —» 3/2~), 8 B(1 + —> 2 + ), 
it is found that their inclusion leads to agreement with the experimental data. This scheme has been 
most recently utilized to study e.m. (both E2 and Ml) transitions occurring in 8 Be [43, 44]. It is 
found that the agreement between the calculated and the experimental Ml widths is not satisfactory. 
Nevertheless, chiral two-body e.m. currents provide correction at the 20%-30% level, which, in most 
but one case, improves the IA values. It is possible that the systematic underprediction of these ob¬ 
servables is due to a poor knowledge of the small components entering the calculated GFMC nuclear 
wave functions 0. 
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5 Summary and outlook 


In this talk, I presented an overview on the present status of ab initio calculations of e.m. observables, 
including e.m. moments and form factors, as well as e.m. transitions in light nuclei. The emphasis 
was on calculations that account for many-body effects in both nuclear Hamiltonians utilized to gen¬ 
erate the wave functions, and e.m. current operators. I focused on results that account for two-body 
operators that have been derived from a ^EFT formulation with pions and nucleons, including up to 
corrections of two-pion range. The ab initio prescription is extremely successful in explaining the 
experimental data, provided that many-body effects in both the e.m. currents and nuclear Hamiltoni¬ 
ans are accounted for. ^EFT based calculations of ,4=2 and 3 nuclei e.m. form factors 113ill nicely 
agree with the experimental data in the low-energy regime of applicability of^EFTs, with two-body 
corrections playing an important role in improving the agreement between the calculated and the ex¬ 
perimental values of the trinucleon magnetic form factors. E.m. two-body current operators provide a 
40% correction in the calculated ma gne tic moment of 9 C |40i], and corrections at the 20%-30% level 
in Ml transitions occurring in 8 Be [44]. There are many interesting e.m. observables that can be 
accessed within this formalism. For example, few (or no) ab initio calculations of e.m. (charge and 
magnetic) form factors in A > 4 nuclei currently exist |j]], and it would be interesting to perform them 
to have a deeper insight on nuclear e.m. structures. A complete microscopic profile of nuclei includes 
also studies of e.m. reactions such as radiative captures and photonuclear reactions. From the theo¬ 
retical point of view, i) the construction of chiral potentials compatible with Quantum Monte Carlo 
computational calculations [45], opens up the possibility of performing consistent Quantum Monte 
Carlo calculations that use chiral potentials and chiral currents; if) the construction of the chiral NN 
potential with the explicit inclusion of A-excitation |46j] , allows study of the effects of A-isobars in chi¬ 
ral two-body e.m. current operators |28|]. Weak processes are also being vigorously studied within the 
^-EFT formulation. Among these studies are the derivation of the axial two-body current operator up 
to one-loop 0, as well as the construction of two-body operators entering pion-production reactions 
induced by neutrino scattering off nuclei (see contributions by F. Myhrer in this proceedings). 
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